MacFormat 1994 August/August CD.bin/Shareware/Education/numericalmethods Folder/chap_7/gauraw.m

Text File | 1994-06-05 | 505 b | 33 lines | [MATF/MATL]

function [A,W] = gauraw
load gauaw;
A = zeros(17,48);
W = zeros(17,48);
c0 = 1;
for j = 1:17,
  m = j;
  i = j;
  if (j > 10),
    m = 12 + 4.*(j-11);
    i = fix(m);
  end
  if (j > 14),
    m = 24 + 8.*(j-14);
    i = fix(m);
  end
  for k = 1:fix((m+1)./2),
    A(j,k) = gauaw(c0);
    c0 = c0+1;
  end
  for k = 1:fix(m./2),
    A(j,i+1-k) = - A(j,k);
  end
  for k = 1:fix((m+1)./2),
    W(j,k) = gauaw(c0);
    c0 = c0+1;
  end
  for k = 1:fix(m./2),
    W(j,i+1-k) = W(j,k);
  end
end
clear gauaw;